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1. Introduction 

The relativistic Hermite polynomials (RHP) were introduced in 1991 by Aldaya et al. |I] in a generalization 
of the theory of the quantum harmonic oscillator to the relativistic context. These polynomials were later 
related to the more classical Gegenbauer (or ultraspherical) polynomials in a study by Nagel [2]. Thus some of 
their properties can be deduced from the properties of the well-known Gegenbauer polynomials, as underlined 
"0"V ky M. Ismail in [3]- In this report we give new proofs of already known results but also new results about 
these polynomials. We use essentially three basic tools: the representation of polynomials as moments, the 
■ subordination tool and Nagel's identity. 
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2. Definitions and Tools 

2.1. Polynomials. 

2.1.1. Relativistic Hermite polynomials. 



. The relativistic Hermite polynomial of degree n and parameter N ^ is defined by the Rodrigues formula 

/ v2 \ N + n Jn / v2\ ~ N 

arm _(-!>•(!+£) 

examples of RHP polynomials are 

N I v\ _ 1 . TjN I v \ _ n v . TjN/ v \_n 1 . v 2 



^ I H§> (X) = 1;H? (X) = 2X; H? (X) = 2 -1 + X' 2 



1 

TV 



^(X)= 4 (l + i) (^( 2+ l)-3X 

These polynomials are extensions of the classical Hermite polynomials H n (X) that are defined by the 
■ Rodrigues formula 

ff„(X) = (-l)%xp(X 2 )^exp(-X 2 ) 
I/"") ' and thus can be obtained as the limit case 

^h"; lim H^(X)=H n {X). 

An explicit formula for the relativistic Hermite polynomial is [1] 

>: m , h n (x] (2N) n ^ (~l) fc n! 

h ; 

where (n) k is the Pochhammer symbol. 
2.1.2. Gegenbauer polynomials. 

The Gegenbauer polynomial of degree n and parameter N is defined by the Rodrigues formula 

C (*) = ^ ("I)" (1 - * 2 ) ^ (1 " A 2 )"^ 

with 

..N (2N)n 



In 



2 n n\ (iV+i) 



examples of Gegenbauer polynomials are 

Cf(X) = l; C(X) = 2iVX; C (A) = 2N (N + 1) X 2 - N; Cf (A) = 21V (iV + 1) ( li^L±2lx : '' X 
and an explicit formula is 

L^J (ah 

(2.2) C n " (A) = £ (-1)* (n l _^ fc , (2A)- 2fc . 
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2.2. tools. 

2.2.1. Nagel's identity. Nagel's identity [2j is a first link between relativistic Hermite polynomials and Gegen- 
bauer polynomials. The RHP with degree n and parameter TV is linked to the Gegenbauer polynomial with 
same degree and same parameter via Nagel's identity 



As a consequence of Nagel's identity, we deduce the following theorem 

Theorem 1. The relativistic Hermite polynomials is related to the Gegenbauer polynomial as 



(2.3) Ci? (X) = a^HZ- N - n -iX X \ -N-n 



2 



with 



a 



N 



(-2*)" 



(2N + n) 

Proof. By the identity (7.2.15.8)], 



•'•V , ; v > _ •>•>" (^)n /1 i \-2\ f rih- N ~ n I X 



c% (ix) = {-2i) n t^t^V i + x 2 ) 2 cl 



(2N + n) n ^ ' VvT+X* 

and the result follows by application of Nagel's identity. □ 

Formula (12. 3p was derived by Nagel [2], who notes: "This representation does not seem to be very useful how- 
ever." Indeed, the fact that the parameter TV of the Gegenbauer polynomial is transformed, for the Relativistic 
Hermite polynomial, into a parameter \ — N — n that depends on the degree n, seems to make this identity 
a priori less useful than expected. However, as will be shown in Section this formula allows in fact an easy 
extension of the scaling identity for Gegenbauer polynomials to the case of Relativistic Hermite polynomials. 

A similar version of formula l|2.3p appears in [8] and also in [9] in the case of deformed Hermite polynomials. 

2.2.2. The subordination tool. Writing polynomials as scale mixtures of others allows to deduce properties 
between them. We use the probabilistic expectation to denote the subordination by the mesure / as 



E b H n (xVPj = £ H n [xVb 



f(b)db. 

The subordination dependence between Hermite, relativistic Hermite and Gegenbauer polynomials is as follows. 
Theorem 2. The Hermite, relativistic Hermite and Gegenbauer polynomials are related as follows 

N% „ „ N ( XVn\ „ , n + 1 



(2.4) C? (X) = -^E b H n ( xVbj , b~T 



n 



N , 
2 



Here, b ~ T [N + §) means that b is a random variable that follows the Gamma distribution with shape 
parameter N + § that is ~/ N+ ™ (6) = w 1 „, e~ b b N+ ^- 1 , b>0. 

Proof. Since 

n+V 



Ec = N 



we deduce from (12.11) that 



,,XVN\ (2N) n LjJ ( _i)fc „| / n+1 \ 

**»( — ] - (^r + i) fc (»-*)!« (^ + — J. ~ (2X) 



(7V + £ 



■An PO 



(2VN) 

so that the first result follows after application of Euler's duplication formula. 
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The same way, we compute 



(N)„ 



2 -E b H n ( xVb 



LfJ 

V (-i) fe 

k=0 y ' 

^2 (-i) fe , (JV) "~*, , (2x) n - 2k 



y(^ + f)^ fe (2^r 2fe 



fc=0 



(n - 2k)\k\ 

which coincides with (12. 2p . □ 
2.2.3. The moment representation. The well-known moment representation of the Hermite polynomial 
(2.5) H n (X) = 2 n E z (X + iZf , 

where Z is Gaussian centered with variance \ is a consequence of the integral formula 8.951]. 

Its extension to the Gegenbauer polynomials is deduced from integral [5j 8.931] (called Laplace first integral 
in |3J) and reads 



(2.6) 



X + - X 2 Z_ 



where the random variable Zn has a Student-r distribution 



(2.7) 



fz N (Z) 



N-l 



1< Z < +1. 



r(A0r(i) 

Using Nagel's identity, we deduce from (|2 .6(1 the moment representation for RHP polynomials 



(2.8) 



X 



iZ 



N 



with Zn distributed according to l|2.7p . 

Another moment representation for the Gegenbauer polynomials is as follows: 

Theorem 3. The Gegenbauer polynomial and relativistic Hermite polynomials have moment representation 



(2.9) 

and 

(2.10) 



2 n (N)i 



-E. 



t,b (XVb- 



iZ 



H 



N 



xVn) = 



2" (2JV) n 



N2 



E bz (xVb + iVT+X^Z 



where Z is Gaussian centered with variance | and independent of b which is Gamma distributed with shape 
parameter N + t-. 

Proof. This expression is derived by the application of the subordination identity 112.41) to the moment repre- 
sentation (12.51) . The representation (12.10p is deduced from (12.91) using Nagel's identity. □ 



However, another set of moment representations involving two random variables can be deduced from (12. 6p 
as follows: 



Theorem 4. The Gegenbauer polynomial has for moment representation 



(2.11) 



{X) = —^e [(x + Vx 2 - l) u + (x - Vx 2 - l) 



-IV 



where U and V are independently distributed according to a Gamma law with shape parameter N. 

Proof. Consider U and V independently distributed according to a Gamma law with shape parameter N; then 
U + V is Gamma distributed with shape parameter 2N so that E (U + V) n = (2N) . We deduce from (|2.6p 
that 



C% (X) 



E{u + vy 



Ep 



X 



i\J\~X 2 Z_ 



N 



but a well-known stochatic representation for Zn is 
(2.12) Z N = 



n 


1 




n\ 


u - 


V 


u + 


V 



X (U + V) + i\J\ -X 2 Z N (U + V) 



where Zn is independent of (U + V) so that 

1 

and the result follows 



Gn{X) = - ] E X{U + V)+Wl- X 2 {U -V) 

TV. 



□ 
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From this result we deduce 
Theorem 5. A moment representation for the relativistic Hermite polynomial is 

(2.13) H% (xVN) = —^E [(< + A) U + (-* + X) V] n 

where U and V are independently distributed according to a Gamma law with shape parameter N. 

Proof. This is a direct consequence of Nagel's formula. □ 

Remark 6. Representation (I2.8[) can proved directly from the explicit expression l|2.ip of the Relativistic Hermite 
polynomials: since Zn has odd moments equal to zero, we have 



the even moments can be computed as 



EZf = 



r( fc + i)r(iv + i) _ 2k! 



T(N + k + l)T(l) klV»(N+$) k 
so that 

-n^ EzAX+iZn) = jf^yh^+h) k ^-2k W [tw) 

which coincides with (12. ip . 

We note moreover that the moment representation (12. lip was derived recently by Sun 0] using another proof 
based on the generating function of the Gegenbauer polynomials. 

2.2.4. example of application. As an example of the usefulness of the moment representations given above, we 
derive the famous [5j 8.952.1] 

4=H n (X) = -^-2 n E (X + iZf = n2 n E (X + iZ)^ 1 = 2nH n _ 1 (X) , 
dX dX 



and its relativistic version 



d rN(x) d (2N) fX 



(2iV) « n F ( X +i7 



N n/2 ^ '* 



n(2N + n-l) IjN 



N 

In the Gegenbauer case, we rather use the stochastic representation l|2.4p to obtain 

ix c " w = ^r^+t jx ( H - ( xVi )) = ^r 2 ^- N+t vhH^ (xvi) 

and since E b ^r N+ ^ Vbf (6) = ^^j^^ E c~r NA „+i / (c) , we deduce 

* C » (X) = ^2/(^+1 + ^ c^V (X) = 2NC^ (X) 

which coincides with [5l 8.935.2]. 

3. The Gram-Schmidt operator 

A family of orthogonal polynomials can be obtained by applying the Gram-Schmidt operator to the canonical 
basis {1, A, ... , A™} . We show here how this operator can be expressed in the case where a moment formula 
exists. 

Theorem 7. If a polynomial P n (A) can be expressed as 

P n (A) = E [X + iZ] n 

for some random variable Z then 

p »w = 

where 4>z (u) = Ez exp (iuZ) is the characteristic function of Z. 
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Proof. By definition 



+ °° -k Jk 

E[X + iZ] n = ( " ) i k EZ k X n - k = TjEZ k ^;x 



k=0 



k=0 



dX k ' 



E z exp[iZ—)X> 



□ 



As an application of this theorem, we recover the following well-known result for the Hermite polynomials 

(2X) n . 



H n (X) = exp - 



4dX 2 

The extension of this result to the case of the Relativistic Hermite polynomials is as follows 

Theorem 8. The Gram- Schmidt operator associated to the relativistic Hermite polynomial is 

d 



H» (XVN) 
where the normalized Bessel function is 

(3-1) j N+i (u) = 2 N +^r(N + 



5 \ dX 



x n 



3\ J N+i («) (-if 



N+- 



2 A- 



* -o kl ( N +Dk 

Proof. The characteristic function of the random variable Zn in (|2 .8[) is the normalized Bessel function 

<S>z N (u) = j N+ i (u) 



□ 



As N — > +oo, since 



(2A0, 



lim 

jV->+oo N~ 



■j N+ i (u) = exp - 



we recover the classical Hermite case. 



4. ADDITION THEOREMS 

4.1. A new proof of the classical addition theorem for Hermite polynomials. The summation theorem 
for Hermite polynomials [5, 8.958.1] states that 



(4.1) 



,2\ ~ / 



-Hn 



| J2k=l a k X k 



e n ™ 

^— ' J - J - m^! 



mH yra T —n k—1 

We give here a short proof using the moment representation; we assume first that J2k=i a k — 1 so that we 
expand 



2 2 "£ 



Zi,...,Z r 



E «fe (-Xfc + iZ k 



\k=l / Lfc=l J lk=l 

where variables Z k are independent and Gaussian with variance | so that X^=i a kZk is Gaussian with variance 
\ and we deduce 



\fc=l / mH hm r - n A;— 1 



r m fc0 _2 



-H mk (X k ) 



so that the result follows. Now we replace a k by , ° fc in this equality so that we obtain the general case 

631. 

4.2. An addition theorem for the relativistic Hermite polynomial. From the well-known addition for- 
mula [7, 7.2.13.36] for Gegenbauer polynomials 

n (N) 

C% (X + Y) = E 7 _7t {2X) n - k Cf + "- fe (Y) 



we deduce the following 

Theorem 9. An addition theorem for the re 



Hermite 



s is 



k=Q 



H» (X + Y) = J2( n ) (1-2N- n) n _ k H? (Y) 
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Proof. Using formula (|2 .3[) . we deduce 



a n m {-iX - lY) = ^ ^ _ j^j (2X) a k H k 2 {-iY) . 



Replacing X by iX and Y by iX and computing 



<*? + "~ fc ( _ 9i sk-n nl(N + n- k) k (2N + n) n (§ - N - n) ' 

„.n — \ Li ) 



k!(N) n (2N + 2n-2) k 

r o^-r* g| r(27V + 2n-fc)r(jV) /l _ 
1 l) k\T(N + n-k)T(2N + n)\2 



so that, replacing X by and Y by iY yields 

Ht n - N (x + Y) = ±(^ (- x | (2^ + <_ fc i?r Ar -"(y) 



2 

and the result is obtained by replacing ^ — TV — n by TV. □ 



5. The scaling identity 
The scaling identity for Hermite polynomials reads [3J 4.6.33] 

LtJ , 

(5-1) H n (cX) = £ (-1)' - (1 - c 2 ) 1 c»- 2 ^„_ 2i (X) 



A quick proof can be given using the moment representation (|2,5p : 

H n (cX) = 2"£ z (cX + iZ)" 



= 2"^,^ [cX + icZx + iy/l - c 2 Z 2 
where Z\ and Z 2 are independent Gaussian random variables with variance h so that 

H n (cX) = 2 n J2 (l)i k (l~c 2 y EZ*c n - k E(X + iZ 1 ) n - k 



k=Q 



LfJ f \ 

= 2 " E (2/) ^ ( X - c 2 ) 1 EZ 2 2 l c n - 2l 2 2l - n H n _ 2l (X) 

since the odd moments of a Gaussian are null; as moreover EZ 21 = 2 2l ~ 1 r ^ 2 ^ we deduce the result. 

It is possible to extend this proof to the case of Gegenbauer polynomials using either the moment repre- 
sentation l|2.6p or (I2.9[) : however, a more simple proof can be derived using the subordination relation (12.411 as 
follows. 

Theorem 10. 7.2.13.37] The scaling identity for Gegenbauer polynomials reads 

(5.2) C» (aX) = £ tllMl (1 _ c2) < c^C^ (X) 

1=0 

Proof. From the subordination formula l|2.4p and the scaling formula (15.111 we obtain (where the notation &at + « 
is a shortcut for b ~ r^a) 

1 LfJ , 

E b H n (cX^b~^) = ^-C» (cX) = £ (-1)' 7-^)77! I 1 " ^E b U n . % [xjb^ 

but 

£# n _ 2i (X^/^T) = EH n , (x^b N , + ^ = jj^C$ (X) 
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with ri = n - 21 and N' = N + I so that 



C (cX) 



(N)„ 



LfJ 



n 

LfJ 



(rc-2Z)!Z! 



(1 - c 2 ) 



' n -2i (n-2l)\ N+l 



C„-2i (*) 



The scaling identity for the RHP can be deduced from the preceding one using formula (|2,3p . 
Theorem 1 1 . The scaling identity for relativistic Hermite polynomials is 

LfJ , 



□ 



1=0 



Proof. Starting from (|5,2p and using l|2.3p . we deduce 



i-N-T 



,N+l 



-i(X\l- -N-n 



L f J ( W 



;=o 



2\ ; „n-2i 77-5^^^™+' 
C -"n-2Z 



-/.A'W- -iV-n + Z 



A short computation yields to 

^n-2; 



(|-iV-n + Z) 2 T(iV + n+i) 



1 k < 1 ^ (n-2Z)! (I_7V_ n )f T(N + n-l + l)- 



Replacing TV by | — n — N we deduce 

(cXVN) = ±^ l)m N , r(1 _ N _ t) 
and the result follows from the fact that 

r(1 "^ -(-i)' m 



£3 {xVn 



+ / 



r(i-jv-z) 



□ 



6. Generating Functions 



6.1. the generating function for the RHP. The generating function for the RHP is computed in [6| using 
a differential equation; we note that it can not be obtained directly using formula l|2.3p . However, it can be 
easily obtained from the moment representation i|2.13p or (|2 .8(1 as follows. 

Theorem 12. the generating function for the RHP reads 



E 



/" = l ll-Sl 



N 



-N 



N 



for \t\ < 



/N 



1+- 



Proof. Starting from (|2 . 13[) we obtain 



H n ( X ) t n = E WE) 



n=0 



n=0 



i + 



X 



X 



u - 



U Ey exp 



X 



V 



Ejj exp 



with Euexp (AC/) = (1 - A) w for |A| < 1 so that, for \t\ < , 



X 



+00 N , , 
n .! 



n=0 











tx\ 2 


4) 







-N 



-i + 



X 



V 



-N 



-N 
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We remark that we recover the generating funtion for Hermite polynomials as N — > +00. The proof using the 
moment representation (|2,8p reads 



-2N 



n=0 



1 ( X M V 
7= — F= + lZ N 



t ( X 
= —7= + iZ N 

N \JN 



for \t\ < , But from M 3.665.11 

tp I ury \-2N ^(iV+i) f n 2N . 2 N-l , 2 . ,2\- w 

E ZN {a-ibZ N ) = FTTTwtn / (a-*&cosx) sin zcfe = (a +6) 
so that the result follows with a = 1 — ^ and 6 = -4= . □ 

We note that formula [SJ 3.665.1] can be reovered using only probabilistic tools as follows: since (1 — u)~ 2N — 
Ew 2N ex P (uW) is the moment generating function of a Gamma random variable W2N with shape parameter 
2N, we deduce 

E Zn (1 - bZ N y 2N = E Zn ,w 2N ex P {bW 2 NZ N ) . 

But, by (EH , 

W 2 jv^jv = tr - v 

where [/ and V are two independent Gamma random variables with shape parameter equal to N. Thus 

Ew 2N ,z N exp {bW 2N Z N ) = E uy exp (bU) exp (-6F) = (l - 6 2 ) _JV 

and the result follows. 

6.2. Feldheim and Villenkin. We give here a short proof of the Feldheim Villenkin generating function for 
the normalized Gegenbauer polynomials 

^ C N (1) nf = 6XP (r C ° S ^ jN -i (r sm ^ 

n— n 

where the function ijv-i is defined as in (|3.ip . by remarking that, using the moment representation (|2.6p . 
so that 

51 /-I s ! f = ^ Zjv ( cos0 + iZ N sm9) n — = exp (r cos6>) £ Zjv exp (irZjy sin6>) 



C^(l) n! > n \ 

" v 7 n=0 

The latest expectation is nothing but the characteristic function (\>z N (u) = jjv-| ( u ) of Zjy computed at 
u = rsha.6, so that the result follows. 

An equivalent result for the RHP is as follows: 

Theorem 13. A generating {unction for the relativistic Hermite polynomials is 

+ °° N% „„ / r" 



/V 2 / \ r* 



n=Q y ' n 

Proof. The proof follows the same line as the one above, starting from the moment representation (|2.8p . □ 



* — » rj.l 



6.3. Another generating function for the Hermite polynomials. A classical generating function for the 
Hermite polynomials [3, 4.6.29] reads 

H n +k 



n=0 

where <j> (A, t) — exp (2Xt — t 2 ) is the generating function of the Hermite polynomials. 
A generalization of this formula to the relativistic Hermite polynomials reads as follows. 

Theorem 14. For the relativistic Hermite polynomials, 

£ ^r 1 * 71 - (*. (x (1 + t 

n=0 ' V V / 

where 4>n {X, t) = ( 1 — + X N \ + J is the generating function of the relativistic Hermite polynomials. 



OLD AND NEW RESULTS ABOUT RELATIVISTIC HERMITE POLYNOMIALS 



9 



Proof. Denote 



so that 



/(*,*) = 1-2 



Xt X 2 t 2 



N 



t z 



N 2 N 2 



,i\ 

n=0 

and, using the moment representation (|2 . 13[) . 
^ H» +k (X) 1 ft 



+ OO JV - „s 



Eu,v 
d k 



i + 



X 



U 



X 
N 



i k 



V 



exp 



Iv 
t 



V 



n-\-k 



x 



u 



^jv {X,t) 



Define /3 = yjl + ^ and z = 1 (t - ^) so that (j) N {X, t) = p 2N (l + ^ 

2 x -N 



-N 



dt k 



4> N (x,t) = p 



2 A 



d k 



1 



N 



— p 2N + 2k (— 1) K I 1 



and 

„2\ - N - k 



N 



X 
N 



V 



N+k 



1 + ^ 

' N 



A? 



-JV-fc rrjV 



so that the result holds. 



□ 



7. Determinants 

Determinants with orthogonal polynomials entries have been extensively studied [11] by Karlin and Szego 
and recently revisited by Ismail [12]. We show here that, in the case of Turan determinants, the moment 
representation derived above allows to extend some of these results to relativistic Hermite polynomials. 

We propose a method slightly different from the one used in [12] based on the following result. 

Theorem 15. // the polynomials P n (X) can be expressed as 

P n (X) = E Z [X + iZ] n 
for some random variable Z then the Turan determinant 

Po(X) 

D p n (X) = det 



Pn(X) 



Pn{X) 
Pn+1 (X) 

P2n (X) 



is a constant equal to 



(n + 1)! 



-Ez ,...,z, 



n Vi-z*? 

0<j<k<n 



Proof. We use the formula of Wilks [13j 



det 



to 
mi 



m n 
m n+ i 

m 2 „ 



(n + l)\ Eu °'-' Un II M ~ Uk 



0<j<k<n 



where mu — EUq and the Ui are independent and identically distributed. Thus since 



P n (X) = E z U n 



with U = X + iZ, we deduce 



Uj -U k = (X + iZj) - (X + iZ k ) = 1 {Zj - Z k ) 
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so that 

D « ix)= 7n+i'v n ^-^) 2 

and the result follows. □ 

Using this result, we deduce the Turan determinant for the normalized Relativistic Hermite polynomial 
defined as 

(7.1) K{X) = ^yH? l {x^N 

Theorem 16. The Turan determinant for the normalized polynomial 17.1)) is a constant equal to 

■jCn+l) 



U n y A >- 9 „(„+i) 11 



2 „(n + l) iX^.^.^+l). 
Proof. From Theorem [15] we have 



V 0<i<fe<n 

where, from (|2,8p . each 2j is distributed according to l|2.7|l . This expectation is a Selberg integral equal to |14l 
17.6.2] 

and the result follows after some elementary algebra. □ 

We note that this result can be proved as a consequence of theorem 5 by Ismail [12]: for normalized Gegen- 
bauer polynomials defined as 



C " W= (2A0" C " (X) ' 



the Turan determinant equals 



Applying Nagel's identity, 



so that 



(X) = (1 + Xr (n+1) Di (X) (^j=^) ■ 
Elementary algebra yields the result. We remark the similarity between the former formula and Nagel identity. 

8. Conclusion 

Some new results about Relativistic Hermite polynomials have been shown; the important fact is that several 
tools (subordination, moment representation) have been used, depending on the type of the result. 
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